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If a Z' gauge boson of a gauged L M — L T symmetry is very light, it is associated with 
a long-range leptonic force. In this case the particles in the Sun create via mixing of 
Z' with the Standard Model Z a flavor-dependent potential for muon neutrinos in 
j — . terrestrial long-baseline experiments. The potential changes sign for anti-neutrinos 

and hence can lead to apparent differences in neutrino and anti-neutrino oscillations 
without introducing CP or CPT violation. This can for instance explain the recently 
found discrepancy in the survival probabilities of muon neutrinos and anti-neutrinos 
in the MINOS experiment. We obtain the associated parameters of gauged L^ — L T 
required to explain this anomaly. The consequences for future long-baseline exper- 
iments and for the anomalous magnetic moment of the muon are discussed. The 
main feature of our explanation is that atmospheric neutrino mixing has to be non- 
maximal. Neutrino masses tend to be quasi-degenerate. 
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1 Introduction 

Additional gauged U(l) symmetries are a feature of many theories beyond the Standard 
Model (for a review, see e.g. Ref. pQ). A large amount of interesting phenomenology arises 
in such scenarios, including LHC physics, lepton flavor violation, dark matter, etc. Here we 
focus on a particularly interesting class of models, namely anomaly free U(l) symmetries 
under which the SM is invariant. It was observed long ago [2] that with the particle content 
of the Standard Model one can gauge one of the lepton flavor combinations L e — L^, L e — L T 
or L^ — L T without introducing anomalies. If the gauge bosons associated with this U(l) 
symmetry are very light, then long-range forces are introduced. In case the extra U(l) 
corresponds to L e — L^ or L e — L T , the electrons in the Sun or the Earth generate a 
potential acting on the neutrinos in terrestrial experiments [3H5]. The flavor dependence 
of L e — L M or L e — L T induces modifications to the neutrino oscillations and therefore the 
coupling of the U(l) can be constrained. The lack of a significant amount of muons in 
the Sun or Earth lead to the fact that the oscillation phenomenology of gauged L^ — L T 
with very light Z' was never studied, though this symmetry was analyzed with different 
phenomenology in mind [STTTO]. 

In the present letter we note that the unavoidable Z-Z' mixing in models with gauged U(l) 
symmetries allows to put limits on the parameters associated with L^ — L T . The flavor 
dependent potential generated by the Z' has different sign for neutrinos and anti-neutrinos 
and can therefore lead to seemingly different neutrino and anti-neutrino parameters. We 
apply this to the recently found discrepancy in the survival probabilities of muon neutrinos 
and anti- neutrinos by the MINOS collaboration [TTJ. In this long-baseline experiment, the 
results for the oscillation parameters in the neutrino and anti-neutrino running lead to 
different values, namelju 

Am 2 = (2.35±8;J|) x 1(T 3 eV 2 , sin 2 29 > 0.91 , 
A^2 = (3.36±°;*j) x 1(T 3 eV 2 , sin 2 26 = 0.86 ± 0.11 , ( " ' 

for neutrinos and anti- neutrinos, respectively [TTJ . We will use here the impact of a long- 
range force associated with the Z' of gauged L M — L T to explain this anomaly. We obtain 
the parameters (Z-Z' mixing and gauge coupling) of the U(l) and discuss in addition 
consequences for future long-baseline neutrino oscillation experiments and the anomalous 
magnetic moment of the muon. An interesting feature of our proposal is that in order for 
gauged L M — L T to be the explanation of the MINOS results, atmospheric neutrino mixing 
needs to be non-maximal. We furthermore find an interesting correlation in what regards 
the sign of the differences between neutrino and anti-neutrino parameters. Neutrino masses 
tend to be quasi-degenerate. 

Previous possible explanations for the MINOS results are CPT violation [T2], sterile neu- 
trinos plus a gauged B — L symmetry with a massive (~ eV scale) Z' [13], or non-standard 
interactions [13]. The first two papers [12] and [13] were motivated by previous low statistics 
results from MINOS, while Ref. [13] and the present work use the recent higher statistics 



1 This result is henceforth referred to as "MINOS anomaly" 



data sets ITT]. 



In Section |2] we outline the framework of gauged L^ — L T symmetry including Z-Z' mix- 
ing, current constraints are described in Section [3] The results are applied to oscillation 
phenomenology and the MINOS results in Section HI where we also study the impact on 
future neutrino oscillation experiments, the anomalous magnetic moment of the muon and 
neutrino masses. Section [5] summarizes our findings. 

2 Gauged L^ — L r Symmetry 

The most general Lagrangian after breaking the 577(3) x 577(2) x U(l)y x U(1) L -l t 

symmetry can be written as 



£ — £sm + ^ z ' "*" ^ mix ' (2) 

where the relevant part of the Standard Model Lagrangian is 

£sm = -\b, u b^ - \w; v w^ + \m\ z, z» - ^-f B b, - 4-j% w; , (3) 

4 4 L C\y S\v 

and the hats denote that we are not in the mass eigenbasis. The currents jg and j% are 
the usual Standard Model ones. The gauge coupling of the ^7(1)l m -l t is denoted g'. The 
Z' part in our case is 



C z , = —Z> Z"» + -M' z 2 Z> Z'» - g>r Z' , (4) 



/^nV + ^T^L^-H^-^/^^, (5) 

with the projection operator Pl = |(1 — 75). The term \M'z Z' Z' M breaks the Z7(l)i M _i, T 
symmetry, and is generated by a vev of some Higgs sector (left unspecified here). Then 
there are terms associated with mixing of the field strength tensors and the two massive 
bosons: 

A™ = ~^ Z'^ B, v + 5M 2 z; Z» (6) 

with the kinetic mixing angle x- The crucial mixing term sin x can arise directly, or can 
be generated radiatively [T6] . 

Diagonalizing [15] the kinetic terms (which gives fields denoted by B^ = B^ + sin x Z and 
Z'^ = cos x Z) and then the mass terms leads, besides the usual W bosons, to a massless 
photon field A^ = dw B^ + sw W 3fl and two massive gauge bosons Z\ and Z 2 . They are 
related to the original Z and Z as 

Z% = cos £ [Z^ - s w sin x Z'^J + sin £ cos x Zj, , (7) 

Z% = cos f cos x Z^ - sin £ \Z* - s w sin x Z'A , (8) 



where £ is a new mixing angle defined by 

tan 2f = -, - Z1^XMS1±^MK^4 . (9 ) 

M 2 , - M§ cos 2 x + M 2 s 2 w sin 2 x + 25 M 2 s w sin x 

The above physical particles Z\ and Z 2 are in the literature normally called Z and Z'. We 
will follow this notation from now on. Their masses are given by 



^ '■'-T- ± \l 1 ' 2 (V^ ( 10 ) 



r r> a + c / , „ fa — c 



with 



o = Mi , b = s w tan x Mi + 

1 , . C ° SX . N (") 

c = — — I M§ sL sin 2 x + 2Sh/ sin x 5M 2 + Mf , I . 
cos" 1 x v ' 

The situation simplifies considerably if the Z' is much lighter than the Z, i.e., if x "C 1 
and 5M 2 <C M§ are very small. In this case we have for the masses 

M 2 ~M 2 Z , M 2 ~c- — , (12) 

a — c 

and the mixing angle is 

1 (, . 5M 2 \ „ 5M 2 

£~ \8 W smx + -7r- \~swX + -^r- U 3 

cosx y M\ J M 2 Z 

With this approximation the Lagrangians for the physical particles arq 2 ] 

£a = -e{jjm)nA> J ', 

£zi = -\ ((h)» - s 2 w (j'em) m ) + gftUX ) Z l > 

£z 2 = - ( #' (j")m - (^ - s h/ x) (O's)/* - s w (Jem)/.) - e c w x (Jem)/* J Z£ . 

V s w cw / 



(14) 



The Lagrangian for the A^ 1 field is the canonical one and hence e — e. The other gauge 
coupling g' is simply g'. 

If we take the mass of the Z' to be M 2 < l/R A .v. ^ 10" 18 eV {R A . V . ~ 7.6 x 10 26 GeV" 1 de- 
notes an astronomical unit) we obtain for particles on Earth a static potential generated 



^TToro uro rlofinorl fVio nVn7-cir»cil AA/ninKn-rrr Qnrrlo oq c~__ /^f__ — i LL_ 



''Here we defined the physical Weinberg angle as s^ c w = ,- v *j 2 . This gives the identity 
Sw cvk Mi = sw cw Mz and the neutral current coupling constant becomes e/(sw cw) — e/(sw c w) (1 ~ 

em- 
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Figure 1: Long-range v^ r — (e,p, n) interaction through Z-Z'-mixing. 

by particles in the Sun. This has been studied for the U{1)l £ -l^ and the U{1)l,,-l t gauge 
bosons, for which the electrons in the Sun generate a potential 

^ a ^^- 3 ' irI1 (i^) eV (15) 

for the neutrinos i/p on Earth. Here a e p = g' 2 /{An) is the "fine-structure constant" of 
the U(l)z e -La and N e is the number of electrons in the Sun. The constraints from solar 
neutrino and KamLAND data are a efJ/ < 3.4 x 10" 53 and a eT < 2.5 x 1CT 53 at 3a [3H3]. 
The lack of muons and taus seems to forbid analogous studies of L^ — L T , since its Z' 
does not couple directly to protons, neutrons or electrons. Consequently, to the best of 
our knowledge, there is no limit on a^ T from oscillation experiments. 
However, there is an indirect effect due to the Z—Z' mixing (see Fig. [TJ. For a neutral and 
unpolarized Sun the final result for the potential is (see the Appendix for details) 

e N 

V,,r = ±g'(Z-sw X) 7 j-^~ ■ (16) 

4 s w c w 47tRa.u. 

Looking at Fig. [[J the main features of this potential can be understood as g' and e/(% cw) 
arising from the vertices and (£ — swx) from the Z—Z' mixing (see Eq. fll4p ). The con- 
tributions of the electrons and protons cancel each other, so that finally only the neutrons 
generate the potential. Their total number in the Sun is about N n ~ N e /4 ~ 1.5 x 10 56 . 
The Earth also generates a comparable potential, approximating a static potential at the 
surface, we get 

Karth N n ^ rth Ra.u. 1.8 x 10 51 1.5 x 10 8 nog 

= — ■ ~ — ~ 0.28 . (17) 

Kun AT„ >sun R suvface 1.5 x 10 56 6380 V ; 

Our full potential at the surface of the Earth is therefore: 

V^ T = ±V = ± 3.60 x lO" 14 eV (j^so) with a = g' (^ - s w X ) ■ (18) 

For anti-neutrinos, the sign of V changes. We stress here that the parameter a that we 
have defined is not a "fine-structure constant" as for the L e — L^ or L e — L T potentials, but 



a combination of coupling and mixing parameters. It can in particular be either positive 
or negative. Note further that due to the various factors in V the scale for a = 10~ 50 is 
different than for a e p = 10~ 50 in the cases of gauged L e — L^ or L e — L T in Eq. ( Fl5l) . We will 
use in the following the value give in Eq. ( Tl8l) for a long-range force according to the Earth- 
Sun distance. In order not to completely spoil the successful oscillation phenomenology, 
V should not become too close to Am 2 /E ~ 2.9 x 10~ 12 (GeV/E) eV, where we took for 
Am 2 the mean of the two mass-squared differences from Eq. ([I]). 

The crucial Z-Z' mixing, and consequently the potential f fl6|) . can only be avoided if for 
the Lagrangian in Eq. ([6]) C m [ x = holds, i.e., if both x an d SM 2 vanish. As can be seen 
from Eq. ( fl3|) . a would vanish for SM 2 = 0. In that case, however, one can show that the 
next order term for £ would generate non-zero a ~ g' sw {Mz'/Mz) 2 Xi which is however 
too small for our purposes, as we will see later. In the case x = 0, the mixing angle is 
given by tan 2^ = „~ 2 SM - 9 ■, and a looks as before. 

to J s M 2 -M 2 Z , ' 

If M 2 < 1/-R g ai — lCT 27 eV, with R ga x the distance between the Sun and the core of the 
galaxy (R ga ,\ — 1.6 x 10 9 -Ra.u.)j we would obtain a potential 

l^(l-4)xl0" 

VU 1-6 x 10 9 y ' 

(with 100 — 400 billion stars) which would dominate over the Earth and Sun potentials. 
Depending on the range of the U(l) force the results which we obtain in the following can 
be easily rescaled. 

3 Current bounds on L^ — L T parameters 

In this Section we will discuss the current bounds on the parameters of L^ — L T . They arise 
from gravitational fifth force searches, electroweak precision observables, fermion charge 
universality and cosmological considerations. 

In principle our model violates the equivalence principle because it adds a lepton number 
dependent force to gravitation. The bounds on such forces are very strict [17] but are 
not directly applicable here since they are based on lunar ranging and torsion balance 
experiments, which are only sensitive to the electron and baryon content. The only effect 
comes once again from mixing; as shown in the Appendix, the potential corresponding to 
Z' generated by a massive body depends on its neutron number iV n : 

v{r) = 'JL^X> N .*p.. (20 ) 

4 sw cw 47r r 

The gravitational potential between two bodies with masses mi and m 2 and neutron con- 
tent N ni and N ri2 is therefore changed to 

v sr „ (r ) = - G „!^i fi - (-AzieA) 2 «««. i e -* | . (21) 

\ 4 sw cw J "mi m 2 4it& n 



The 95% CL. limits for a neutron dependent fifth force as a function of its range are given 
in [T7] (see references therein for a description of the experiments), where the effect of new 
light vector or scalar bosons is parameterized as 

\r i \ n m i m 2 (-, . ~N ni N ni /A \ 

V grav {r) = -G N 1 + a e ' , (22) 

/i being a test body mass in units of atomic mass unit u and a = ±g 2 /(47r Gn u 2 ) (the sign 
distinguishes between vector and scalar interaction). Comparison with Eq. (12T|) gives the 
translation into our parameters 

\a\ = l f e ^~ s wX) \ 2 \ = J- r2 ox 

1 l ~47rG N u 2 \ 4s w c w J ' M 2 ' K > 

For Earth-Sun range we take the bound |a| < 10~ n , given in [17], corresponding to 

(£-^x)<5xl(T 24 , (24) 

whereas the limit for an Earth range force is given as \a\ < 5 x 10~ 9 , corresponding to 

(£-^x)<KT 22 . (25) 

These are the strongest constraints on the mixing angles. 

The parameters are however also constrained through precision data from electroweak 
observables. Measurements around the Z-pole examine the mass-eigenstate Z\ with mass 
(see Eqs. ( TT0~][TTj) ) M 2 ~ a(l + b 2 /a 2 ), while measurements on W-bosons give values for 
My/ = Mz cw- Therefore the mixing changes the p-parameter of the Standard Model from 
p = M 2 v /(M 2 c 2 w ) to 

/W = (^-) 2 = PT-^- 2 - P (1 " ^) • (26) 

\MiCwJ l + b z /a z 

The current value [T5] is p = 0.9994 ± 0.0009 which gives £ < 0.025. Stronger limits arise 
by reading off from Eq. ( JT4"j) . the vector/axial couplings of the tauon: 

g v ^2s w ---2 g £, g A -> -- , 27 

2 e 2 

where 2 s 2 ^ — | stems from the SM neutral current j% — s 2 ^ j£ M . The asymmetry parameter 
A T = 2 g v g A / ((g v ) 2 + (g^) 2 ) becomes approximately 

^ -► ^m f 1 + i^l^) ee Al M + AAV0 , (28) 

where A$ M = (1 — 4 s 2 ^)/ (1 — 4 s 2 ^ (1 — 2 s 2 ^)) is the value without any new physics. This 
quantity is measured to be A T = 0.143 ± 0.004 (Ref. [H]), while with the central value 



sin 2 6w(Mz) = 0.23116 one expects A$ M = 0.1499. Since the measured A T and A^ 1 are 
of the same order while a nonzero g' £ shifts them in different directions, we will require 
AA T (g' £) to be within the measured error, i.e. AA T (g' £) < 0.004. This restricts g' £ to 
values 



g'£<3.6 x 10" 



(29) 



This limit is stronger than e.g. from the Z-coupling to v^ or the ratio T(Z — > [i + [i )/T(Z — > 
e + e~), where 



^->'0 = 7 ^%-((^) 2 + (id) 8 ) 



12 s 2 r 2 



(30) 



at tree-level, ignoring lepton masses. 



The mixing also changes the electromagnetic behavior, as can be seen from the La- 
grangian (fT4"|) . slightly rewritten and shown only for negatively charged muons (/i), electrons 
(e) and positrons (e + ): 



C 



z 2 



e ( 1 

sw cw V 4 



/r^/i 



+ 



e c^/ x - (f - s w x) 
ec w x~ (£- s w x) 



s w c w 
e 

SwCw 



' IV 



s 2 
6 VK 



e i_ 7 / ae n 



(31) 



e7/?e } Z^ 



In muonium the coupling between positive muons and electrons is modified because there is 
not only photon exchange, but also photon- Z 7 mixing. In direct analogy to the derivation 
of the neutrino potential given in the Appendix, one finds an effective potential 



V, 



■ '• -rA 1 -^-'*" 



(32) 



where Qp = — (C,—swX)(^/^~ s w)/( s w c\v) — cwX- Hence, the result is an effective change 
of the fine-structure constant in systems involving muons (or tauonsjS On atomic scales 
the factor e~ rM2 can be omitted. By comparing the above potential with the potential for 
positronium we find the ratio of the /i + and positron charge 



Q(e+ 



£d 



9 tQi 



i- 



9 -Q P 

e 



(33) 



3 There is an effect quadratic in Qp due to two mixings in systems like positronium or hydrogen, which 
is however way to small to be observable. 



This ratio has been measured via the muonium hyperfine-structure [23] to be 1 with an 
accuracy of 10~ 7 , corresponding to a limit 

g'(3s wX + (l-4s 2 w )Z) <5xl(T 8 . (34) 

Note that, as it should, there is no effect in case of x — C = 0, i.e., when there is no photon- 
Z' mixing. In case di- muonium (a bound state of \i~ and /z + [24]) would be produced, one 
could test the Z' even in the limit of no mixing. 

Another effect the new light Z' would have is a contribution to the effective number of 
degrees of freedom, potentially threatening for instance the success of Big Bang Nucle- 
osynthesis (BBN). Recent BBN measurements as well as other cosmological probes are 
compatible with about one extra degree of freedom [19] . Let us demand that the Z' does 
not contribute. This means for the case of BBN that it should enter equilibrium after weak 
interactions freeze out (T ~ MeV), and requires to consider the process Z' Z' — > v^ T v^ r , 
whose rate goes as (g' 2 / (in)) 2 T . Comparing this to the Hubble rate H ~ T 2 /M P \ gives 
the requirement g' 2 /(47r) ^ 10~ n [20]. A constraint of similar size has been estimated 
from Supernova 1987a [2T]. An upper limit of g /2 /(4-7r) ^ 10~ 18 can be obtained with the 
process 7/x — > Z' //, going with g' 2 /{An) aT, and demanding that Z' is not in equilibrium 
at T = m M [22] . 

As expected, the largest constraints stem from the equivalence principle and BBN. How- 
ever, the small values of the L M — L T parameters required in order to give observable effects 
in oscillation experiments are compatible with these limits. 

4 MINOS and Beyond 

The potential V in Eq. (IT51) generated by L^ — L T is flavor dependent, acts on the \x-r 
part of the system, and has a different sign for neutrinos and anti-neutrinos. Consequently 
it is a good candidate for an explanation of the MINOS results, which seemingly give 
different mixing parameters in the muon neutrino and anti-neutrino oscillations. In a 2- 
flavor approach, the Schrodinger-like equation for neutrinos is (note that we start in the 
mass basis) 

where Am 2 = m 2 — m 2 , is the atmospheric mass-squared difference and vm = (^2, v?) T are 
the mass eigenstates which are connected to the flavor states ^flavor = (^, ^r) T = U um via 
the matrix 

/cos* sin^y 
\-sinti cos 6 J 




E [GeV] 

Figure 2: Difference between the mass-squared differences of neutrinos and anti-neutrinos 
(choosing initial values of Am 2 = —2.48 x 10 -3 eV 2 and 8 = 1) for different values of a as 
a function of energy. 

Here 9 = 823 is the atmospheric mixing angle. 

The Schrodinger-like equation ( )35|) thus contains the Hamiltonian: 



m 



1 (m 2 + 2EV cos28 2EV sin2# 

m 2 - 2 E V cos 28 



J_ TJ mly 
2E V \ ml v 



2E\ 2 E V sin 28 
As we have indicated, Hy is diagonalized by the rotation matrix 

2 77 sin 2# 



4 



(37) 



tf 



v 



cos (/) sin < 
- sin 6 cos ( 



with tan 20 



1 - 2 77 cos 20 



(38) 



We have introduced 77 

the new mass eigenstates um; 



2EV 
Am 2 



. The new mass eigenvalues m 2 , v and m 2 y are associated to 



l^2,v,^3,vj via 



^Af,V 



c/i 



t4 f/ f 



y V M — Uy U ^flavor 



(39) 



Thus, in the presence of the potential V, the mixing angle between flavor and mass eigen- 



states becomes 8 + and Am 2 changes to Am 
parameters are 



m 



3,V 



m 2 , ,/• The exact results for the 



b 2,V 



sin 2 20i 



sin 2 28 



1 - 4 77 cos 26 1 + 4 77 



2 • 



Am v = Am 2 A/l -47/ 003 26* + 4t? 2 = Am 2 



sin 2 20 



sin 2 28 



(40) 
(41) 



v 



For V = the vacuum results sin 2 2$ and Am 2 are obtained. For anti-neutrinos, the 
potential V and hence 77 changes sign, thereby an apparent difference between the oscillation 
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parameters of neutrinos (Amy, 9) and anti-neutrinos (Amy, 9) could arise. Fig.[2]shows the 

difference between the mass-squared differences of neutrinos and anti-neutrinos (choosing 

an initial value of Am 2 = —2.48 x 1CT 3 eV 2 ) for different values of a as a function of 

energy. 

We note here three important properties following from Eqs. ( 140 



first, the effect goes with 77 cos 29, and therefore it is absent if 9 is maximal. In 
this case the oscillation parameters 9 and Am 2 would be the same for neutrinos and 
anti-neutrinos, but with a common offset compared to their values for V = 0. If 
the long-range force mediated by L M — L T is responsible for the MINOS anomaly, 
then the necessary 9 7^ 7r/4 is a possibility to disentangle it from any other proposed 
explanation 



• the second point is that the corrections to the mixing angle and the mass-squared 
difference are correlated. For positive Am 2 and a the correction for sin 2 29 goes in 
the opposite direction as the correction of the Am 2 . Recalling that MINOS finds 
Am 2 > Am 2 we therefore predict for positive Am 2 and a that sin 2 29 > sin 2 29, 
which is compatible with the MINOS results (see Eq. (JT|)), and can be checked with 
higher statistics data sets. For negative Am 2 and positive a the correction goes in 
the same direction, and hence sin 2 29 < sin 2 29; 

• the third point is that the relative effect is expected to be slightly larger for sin 2 29 
than for Am 2 . 

We can estimate the magnitude of the parameter rj as 

" ^ ° 025 (ws) {Sv) • < 42) 

which allows for not too high energies (note that at MINOS the oscillation dip occurs at 
around E ~ 1 GeV) and for a around 10 -50 (see the discussion after Eq. ffTB"]) ). 1] is small 
and can be used as an expansion parameter. As can be seen from (T40|) and (T4T1) the relative 
difference of the mass-squared differences is in this case obtained as 



Am 2 
while for the mixing angle the result is: 



4r]cos29, (43) 



Am?, — Am v 

1 *__ r^i 



sin 2 29 v - sin 2 29 



y ; bm ^ czSri^e. (44) 

sin 29 

These expressions nicely confirm the three points mentioned above. The muon neutrino 
and anti-neutrino survival probabilities are 

P = P(^ -►!/„) = 1 - sin 2 29 v sin 2 ^L , (45) 

P = P(F„ -> 77^) = P(v^ -> Vli ){a +* -a) , (46) 
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Figure 3: The oscillation probabilities for the best-fit values from Eq. (1491) for neutrinos 
and anti-neutrinos superimposed on the MINOS data. Also plotted are the cases a = 
and the value for the second, local x 2 -minimum. 

which are subject to the following degeneracies 

P(9, Am 2 , a) = P(9, -Am 2 , -a) = P(9 + vr/2, Am 2 , -a) = P(9 + tt/2, -Am 2 , a) . (47) 

While the part discussed so far was rather general, we continue by applying the formalism 
to the recently found MINOS results |llj . We have performed with the expressions ( 1431 
[461) a x 2 -fit to the MINOS data (given in bins of energy Ei) on the ratio of observed events 
divided by the expectation for no oscillations. This data was taken, as in Ref. [H], from 
the slides of the talk referred to in our Ref. pj] . In case of asymmetric errors, the largest 
one was used and inserted in the x 2 -function 



X 2 (#,Am 2 ,a) = ]T( 



P{6,Am 2 ,a,Ei)-R t 



o~ 



+ E 



P(6, Am 2 , a, E { ) - R 



-^ x 2 



(4* 



where P (P) is the survival probability Piv^ — ¥ Vy) from Eq. ( I45j) (from Eq. ( 146]) ). Ri 
(Ri) the ratio of observed events relative to the no-oscillation expectation, and <Ji (Wi) the 
error for the neutrino (anti- neutrino) data set. The result of our fit after marginalizing 
over Am 2 and 9 isj 



sin 2 29 = 0.83 ±0.08 , Am 2 = (-2.48 ±0.19) x lO^eV 2 



a 



(1.52+^) xKT 50 , (49) 



with XmiJNdoi = 47.77/50 ~ 0.96. Recall the degeneracies listed in Eq. (H7j) . In Fig. |3] 
we show the experimental data together with the results of our fit. One can see that the 



4 We have checked our analysis by setting a = and have obtained the best-fit values Am 2 = 2.28 x 10 3 
eV 2 , sin 2 29 = 0.94 for the neutrino data set, and Am 2 = 3.38 x 10~ 3 eV 2 , sin 2 29 = 0.81 for the 
anti-neutrinos, in good agreement with the MINOS results. A fit to the total data set yields Am 2 — 
(2.38t°; 2 ?) x 10-3 eV 2 and sin 2 29 = 0.89±g;^, with Xmin/^dof = 49.43/51 ~ 0.97. 
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Figure 4: The x -function from Eq. ( 1481) as a function of the fit parameters. 



X 2 contours at 1 , 2 and 3cr 



X 2 contours at 1 , 2 and 3<x 




sin 2 26 

Figure 5: The 1,2 and 3<r (or equivalently A% 2 = 2.3,6.18, 11.83) contours in the a-Am 2 
plane when marginalized over 9 (left) and in the a-sin 2 29 plane when marginalized over 
Am 2 (right). The cross marks the best-fit point. 

non-zero value of a puts in particular the data points at the oscillation minimum in better 
agreement with the curves. From the plot of the x 2 -function in Fig. H] one sees that there 
is a second (local) minimum, corresponding to sin 2 29 = 0.98, Am 2 = 2.36 x 10~ 3 eV 2 and 
a = 4.41 x 10" 50 , with Xmin/^dof = 48.73/50 ~ 0.97. The curves for this point are also 
plotted in Fig. [3j The second local minimum also explains the "rabbit head looking" shape 
of the contours in a-Am 2 and a-sin 2 29 space shown in Fig. [5j 

The goodness of fit is not particularly worse for the absence of new physics, which has 
been noted also in Ref. [14] . 



We continue by discussing the consequences of the implied value of a in future neutrino 
oscillation experiments. We have modified the commonly used GLoBES software [22] to 
include the potential V from Eq. ffTBI) . Using the pre-defined packages ("AEDL files") for 
the most frequently discussed future experiments, we analyzed T2K, NOz/A and a neutrino 
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Figure 6: The 1, 2 and 3<r limits which can be obtained by T2K (left), NOz/A (middle) and 
a neutrino factory (right). 
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Figure 7: The 1,2 and 3<r constraints on a which can be obtained by T2K (left), NO^A 
(middle) and a neutrino factory (right) if a = (1. 521^2}) x 1CT 50 . 



factory, as listed in Tabled! to obtain future constraints on a. The oscillation parameters 
we use are listed in Table [2j The result is that at 3a, a can be constrained to be below 
11.80 x 1CT 50 , 1.93 x 10" 50 and 0.53 x 10~ 50 , respectively. The x 2 -functions generated by 
GLoBES are shown in Fig. EJ 

Setting the true parameter values of a, 9 and Am 2 (and their errors) to our best-fit values 
from Eq. ( 149]) . we can see how the "precision" on a can be improved. From the plots of 
X 2 in Fig. Hone sees that NOz/A would give a = (1.52 ± 0.27) x 10~ 50 , T2K would yield 



a 



1.52±0.46) x 10 and NuFact would determine very precisely a 



;i.52j£ji)xio- 



-50 



As mentioned above, long-range forces generated by L e — L^ r have been discussed before. 
Ref. [3] bounds a e/ijT by analyzing v^ and v T oscillations and using atmospheric neutrino 
data. It is easy to see that in a two-flavor framework, the potential V eT = a eT N £ /Ra.v. 
corresponds to 2 V^ T . Likewise, V eil corresponds to —2 V^ T . Therefore, the limit of a eT < 
6.4 x 10~ 52 obtained in Ref. [3] corresponds to a — g' (£ — % \) < 8-9 x 10~ 50 , not in 
conflict with our fit-result from Eq. (H9l) . In turn, this means that not only L M — L r could be 
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the origin of the MINOS anomaly, but also L e — L^ or L e — L T , for which a = 1.52 x 10~ 50 
translates into a efl)T = 1.1 x 10~ 52 . If we take the 3<r-bound a eT < 2.5 x 10 -53 from 
solar neutrino and KamLAND data [4] and treat it like in the 2-flavor case we obtain 
a < 3.5 x 10~ 51 . However, the interplay of the other limits on long-range forces, and also 
the impact of stronger bounds on a eM . T using solar and KamLAND data [3], can not be 
used without doing a full 3-flavor fit to all data. In general, we note that the different 
flavor structures of the potentials arising from L e — L^, L e — L T and L M — L T , 

(50) 

render it difficult to translate existing bounds on L e — L M or L e — L T into constraints on 
L^ — L T , in particular if in addition a matter potential is present in V ee . We would like to 
stress though that the solar neutrino oscillations should really be fitted specifically for this 
model, since the electron and neutron densities in the Sun are not proportional. 
In Ref. [H] the presence of Non-Standard Interactions was assumed as the reason for the 
MINOS anomaly. In particular, a term e Mr was introduced, and in the Hamiltonian it 
appears together with the potential V m — y2Gpn e ~ 1.1 x 1CT 13 eV. By fitting the MI- 
NOS data, the value e MT = —(0.12 ± 0.21) was obtained. We note that for |e MT | = 0.1 the 
term V m e^ T is of the same order of magnitude as our potential for a ~ 10 -50 . A small 
difference to our explanation is that V m oc n e , i.e., the potential depends on the electron 
density, which in turn depends on the matter density of the Earth. This changes with 
baseline, and hence could in principle be used to distinguish Non-Standard Interactions 
from our explanation. We should note here that in a 2-neutrino framework the relation 
2 V = V^e^ holds, and our range of a would correspond to e MAt ^ 0.25, to be compared 
with the 90 % C.L. limit [27] (e^l < 0.068. Saturating this limit would corresponqj to 
a = 1.04 x 10~ 51 . A fit to the data fixing it to this value yields sin 2 29 = 0.88+o;q7 and 
Am 2 = (-2.39+£ 2 ?) x 10~ 3 eV 2 , with xLJ N doi = 49.25/51 ~ 0.97. 

It is worth discussing the anomalous magnetic moment of the muon, where since many years 
a conflict between theory (i.e., its Standard Model calculation) and experiment exists [T8] . 
The current experimental value of a M differs by 3.2a from the Standard Model prediction, 
although there is some uncertainty in the hadronic contributions. Nevertheless, since the Z' 
couples to the muon, it contributes to Aa M [B]. In the limit of M' z <^C m^, the contribution 
is 

A <V = £~ 2 , (51) 

which in our light case translates into a constraint on the coupling g'. From the constraint 
Aa M < 255 x 10' 11 it follows that g' < 4.49 x 10" 4 . This would imply (£-s w x) £ 3.3 x 10~ 47 
in order to explain the MINOS anomaly. 



5 Similar comments apply for limits from atmospheric neutrino oscillations, which have a somewhat 
stronger limit 28 . 
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Turning to neutrino masses, the conservation of L M — L T dictates the effective neutrino 
Majorana mass matrix to be [9|[T0] 



m v = ■ b , (52) 




regardless of its origin, such as some form of see-saw. It would result in neutrino masses 
a and ±b, hence one expects (close to) quasi-degenerate masses. Though the mass matrix 
is fi-r symmetric, and hence implies #13 = and #23 = vr/4, it is too simple and can not 
reproduce all data. Breaking L^ — L T is achieved by introducing extra Higgs particles $', 
which obtain a vev. Necessarily, the implied scale of the Z' mass (which is generated by 
breaking of L M — L T ) and the additional entries in m v are correlated via m' z ~ g' ($') and 
( m ^)a/3 ~ (&') if h is a weak triplet, (m u ) a p & t>„k ($')/A if it is a doublet and couples to 
the SM Higgs, or (m u ) a p ^ ($') 2 /A if it does not. Here A denotes the high energy scale 
which acts as the necessary suppression of the neutrino mass. Simultaneous ultra-light 
Z' of order 10~ 19 eV and sizable (m u ) a/ 3 ~ 0.1 eV implies for, say, A = 10 15 GeV that 
for doublets ($') is of order 10 3 GeV and hence g' ~ 10~ 30 , while g' ~ 10~ 17 for triplets. 
Hence, the Z' will essentially not contribute to the anomalous magnetic moment of the 
muon. 



5 Conclusions 

Long-range forces mediated by the Z' boson associated with gauged L M — L T can lead to 
interesting and largely unexplored phenomenology. For instance, neutrons in the Sun gen- 
erate via Z-Z' mixing a flavor-dependent potential for terrestrial muon and tau neutrinos. 
This potential changes sign for anti-neutrinos, and hence can lead to apparent differences 
in neutrino and anti-neutrino oscillations. Applying this new finding to the recently found 
MINOS anomaly implies a value of around a ~ 10~ 50 , where a = g' (£ — swx) is the 
product of the new gauge coupling and the parameters quantifying the Z-Z' mixing. An 
interesting correlation between the atmospheric neutrino parameters Am 2 and 9 is found. 
The latter is required to be non-maximal, which is one of the handles to probe this ex- 
planation of the anomaly. By making use of the GLoBES software we have furthermore 
discussed future constraints on a. Time will show whether the discrepancy in the MINOS 
results survives. Nevertheless, many new physics effects imply different neutrino and anti- 
neutrino behavior, which underlines the importance of analyzing them separately. The 
new effect arising from L M — L T (via Z-Z' mixing) noted in the present letter is one more 
example for this, and we have given simple estimates for future constraints. 
It would be interesting to discuss a similar approach for other "anomalous" oscillation 
results in which apparent differences of neutrinos and anti-neutrinos are found, such as the 
recent MiniBooNE excess in a z/ M — > v e search [29|, or the slightly larger #12 found in solar 
neutrino analyses with respect to the Qyi m reactor anti-neutrino experiments. As a final 
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remark, neither CP nor CPT violation are in this framework necessary for the apparent 
differences in neutrino and anti-neutrino oscillation probabilities. 

Acknowledgments 

We thank Borut Bajc and Osamu Yasuda for helpful and interesting discussions. This 
work was supported by the ERC under the Starting Grant MANITOP and by the DFG in 
the project RO 2516/4-1 as well as in the Transregio 27. 



Experiment 


Baseline 


Running-time [years] 


Beam-energy [GeV] 


Detector mass 


T2K 
NOz/A 

Nufact 


295 km 
812 km 
3000 km 


5 v +5 77 

3 v + 3 77 

4 z/ + 4 77 


0.2-2 

0.5-3.5 

4-50 


22.5 kt 
15 kt 
50 kt 



Table 1: Parameters of long-baseline oscillation experiments simulated by the GLoBES 

software 



#12 


arcsin V0.318 ± 0.02 (3%) 


#13 


0±0.2 


#23 


arcsin V0.500 ± 0.07 (9%) 


^CP 


G[0,2vr] 


Amlj [10- 5 eV 2 ] 


7.59 ± 0.23 (3%) 


Amli [10- 3 eV 2 ] 


2.40 ±0.12 (5%) 



Table 2: Oscillation parameters [26] used as input to the GLoBES simulation. 



A Derivation of the Potential 



For the sake of completeness, let us give here a derivation of the static potential which 
the particles in the Sun generate for terrestrial neutrinos. The potential (IT5|) for gauged 
L e — L^ or L e — L T can also be derived in this fashion. From Eq. ( 1T4"|) we consider the 
time-like components, note that j^ M = and have that 



■o l-o 1 

J 3 = ~2 6l1 eL+ 2 PL " ; " l>L 



,o 



n L 7° n L 



1 , 

-{n e -n p + n r , 



nn 
' 4 



(Al) 



since the axial-part will result in a spin-operator in the non-relativistic limit and we assume 
the Sun is not polarized. The equation of motion for Z°, following from the Euler-Lagrange 
equation 



<9„ 



S 



Kd v Z 2 ») 



■-Z- 



2a/3 



7 a/3 



5Z. 



2jj 



-MiZ 2a Z« + £ Z2 \ =0 



(A2) 



17 



is therefore 



(«9 2 + Ml) Z° = {Z-s w x)^—^ ■ (A3) 



swcw 4 
In the static case outside of the Sun this is (n n (x) = N n 8^(x)) 



(A - M 2 2 ) Z° 2 = -(£ - SwX )^—- A N n ^\x) (A4) 

sw cw 4 

with the well-known solution 

e 1 e ~rM 2 

V(r) = Z\ = {S-swx) T^n x . (A5) 

swcw 1 * 47rr 

In the limit M2 — > the potential, for v^ and v T respectively, on Earth igj: 

p N 

V„, T = ±g'(Z-s wX h -r^^ + 0{e,x\ix). (A6) 

4 Sw Cw 47TitA.U. 
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